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We present a new method for the quantization of totally constrained systems including general 
relativity. The method consists in constructing discretized theories that have a well defined and 
controlled continuum limit. The discrete theories are constraint-free and can be readily quantized. 
This provides a framework where one can introduce a relational notion of time and that nevertheless 
approximates in a well defined fashion the theory of interest. The method is equivalent to the group 
averaging procedure for many systems where the latter makes sense and provides a generalization 
otherwise. In the continuum limit it can be shown to contain, under certain assumptions, the 
\Q ' "master constraint" of the "Phoenix project". It also provides a correspondence principle with the 

, classical theory that does not require to consider the semiclassical limit. 
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The issue of the dynamics is perhaps the central problem in canonical quantization approaches to totally constrained 
theories like quantum gener al relativity 0,0- There are three salient aspects of the problem that have prevented 
from advancing in the quantization. The first one is how to construct a space of physical states for the theory that 
qq , are annihilated by the quantum constraints and that is endowed with a proper Hilbert space structure. The second 
' issue is related to the introduction of a correspondence principle with the classical theory in particular to check the 
constraint algebra at a quantum level. The third problem is how to address the "problem of time" that is, to 
i introduce a satisfactory picture for the dynamics of the theory in terms of observable quantities. 

We have proposed in previous papers Q a paradigm to deal with the above issues that consists in describing the 
theory in terms of a discrete evolution. This is analogous for instance to lattice QCD, where one uses a discrete theory 
to approximate and define the continuum quantum theory as a suitable limit. In our approach the discretization 
, is carried out in such a way that the dynamics of the discrete theory is unconstrained. We call this approach 
f^i ' "consistent discretizations" . The lack of constraints in the discrete theory bypasses almost automatically the three 
\& , issues mentioned above for the discrete theory, but leaves open the problem of how to define a continuum limit. One 
of the important points of our proposal is that the discretization is carried out in a way that the discretization step is 
dynamically determined and therefore one does not have a direct control on how to take the continuum limit as one 
O" 1 does in ordinary discretizations. 

In this article we want to spell out a constructive technique to define properly the quantum continuum limit. Our 
procedure can therefore be viewed as an alternative to the Dirac quantization procedure for a continuum theory in the 
sense that at the end of the day it yields a quantum continuum theory. We will see that the procedure has attractive 
. i-H , advantages with respect to the Dirac procedure. 

A key element is the introduction of a set of discretizations for a given continuum theory called "uniform dis- 
cretizations" . These are such that the evolution steps are bounded by a value that one chooses in the initial data. 
The value of the constraints of the continuum theory (which are not exactly satisfied in the discrete theory) are also 
bounded throughout the evolution. This is of interest in itself since it is quite non-trivial to find discretizations of 
general relativity for which the constraints remain bounded. The important aspect is that since one controls the 
discretization step through the initial data, one can define properly a continuum limit just by choosing data that has 
a step as small as is desired. 

We start by considering a canonical theory with a phase space with canonical variables qi,Pi that is totally con- 
strained, by this meaning that the total Hamiltonian Ht is a linear combination of N constraints fa which we will 
assume are first class. We are assuming we are dealing with a mechanical system with a finite number of degrees of 
freedom. This is of interest in the context we are discussing since field theories when formulated on the discrete space 
— as is common, for instance, in loop quantum gravity — , become such systems. 

We will now introduce a discrete evolution given by the flux of a Hamiltonian H that is constructed from the 
constraints of the theory in a way we will soon discuss and such that the evolution of any dynamical variable A is 
given by 

A n+1 = e^{A n ) = A n + {A n ,H} + ±{{A n , H}, H} + ■ ■ ■ (1) 

As is obvious, H is a constant of the motion of the discrete evolution. 

The uniform discretizations are given by a family of Hamiltonians H constructed in the following way. Consider a 
smooth function of N variables f(x±, . . . , x n ) such that the following three conditions are satisfied: a) f(xi, ■ ■ . , x n ) = 
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<^> x, = OVi and otherwise / > 0; b) |£(0, . . . , 0) = 0; c) detgg^- ^ Vx and d) /(&(g,p), . . . <f> N (q,p)) is 
denned for all q,p for the complete phase space. Given this we define H(q,p) = f((f>i(q,p), ■ ■ ■ 4>N{q,p))- 

A particularly simple example is H(q,p) — 1/2 ^\ =1 (f>i(q,p) 2 , a choice that has interesting parallels with the 
"master constraint" of the "Phoenix project" Q as we shall discuss later. 

An important point is that if we choose initial data such that H < e then (j>i remain bounded throughout the 
evolution and will tend to zero in the limit e — > 0. Let us see that in this limit one recovers the evolution equations 
given by the total Hamiltonian Ht in the constrained continuum theory. Let H as in the simple example above 
and take its initial value to be Hq = S 2 /2. We define A-j = (pi/5, and therefore = 1- The evolution of the 

dynamical variable q is given by, 

N 

q n+1 = qn + ^{qnAi}h5 + 0{5 2 ) (2) 
i=l 

and if we define q = lim,5_>o (<Zn+i — q n )/8, where we have identified the "time evolution" step with the initial data 
choice for S, one then has, 

N 

q = ^2U,4>i}^i, (3) 

»=i 

and similarly for other dynamical variables. The specific values of the multipliers A^ depend on the initial values of 
the constraints <pi. 

Notice that in the current proposal the evolution is generated by the exponential of a Hamiltonian. In our previous 
treatments of "consistent discretizations" one used as a starting point a type 1 or 2 generating function of a canonical 
transformation. Generically one can associate a generating function with a choice of /, by constructing the Hamilton 
principal function. The knowledge of such generating function could be useful, for instance, to actually use the discrete 
scheme in a concrete numerical implementation classically. As an example, we shall analyze a system with N Abelian 
constraints by using a slightly different technique that seems more convenient for numerical applications. We define 

L(q n , q n +i> -Vii ■ • • : Ajv) = S(q n ,q n +iXi, . . . , \n) + s(Ai, . . . , Xn) (4) 

where L is a type 1 generating function of a canonical transformation between canonical variables q n , p n and q n +i , Pn+i > 
S is Hamilton's principal function for a given set of Lagrange multipliers Ai, . . . , Ajv (they are evaluated at instant n, 

we omit the subscript for simplicity) and g is a smooth function such that g(0) = 0,<9i<?(0) = and det d ° q x ^ 0. 
The generating function yields the canonical momenta in the usual way p n +i = dL/dq n +i, p n — —dL/dq n . One also 
has that dL/dXi = 0. The latter equation determines the Lagrange multipliers, A,; = hi((f>), where hi is the inverse 
function of the mapping defined by Ai — > J^. This evolution corresponds to a Hamiltonian H = f(4>\, ■ ■ ■ , 4>n), with 
dif = hi. In other words, the Legendre mappings induced by g and / are inverse. In particular, when g — y\ —1 xf/2 
then H = YliLi < f ) i/2- The generating function L allows to determine the discrete evolution that preserves exactly 
the value of the constraints of the continuum theory and recovers the continuum limit when all — > in the initial 
data. 

The constants of the motion of the discrete theory are quantities that have vanishing Poisson bracket with the 
Hamiltonian, {Of , H} — and in the continuum limit Hq — > reproduce, as functions of phase space the "perennials" 
of the continuum theory Of — lim# ^o Of . This can be immediately seen from the fact that the discrete equations 
reproduce the continuum equations for any dynamical variable in the continuum limit. Conversely, for every perennial 
of the continuum theory there exists a constant of the motion (in general many constants) of the discrete theory that 
reduce to the given perennial in the continuum limit. We have therefore shown that uniform discretizations recover the 
constraints and the perennials of the continuum theory and therefore provide a good starting point for a quantization 
of the continuum theory. 

We now turn our attention to the quantum theory. We will introduce a Heisenberg quantization for the discrete 
theory (this is more natural given that one has an explicit evolution). To quantize the theory we follow several steps. 
We start with the classical discrete system constructed as in the previous section, we eliminate the canonical variables 
at level n + 1 in terms of the variables at level n, q n +i = qn+i(qn,Pn), Pn+i = Pn+i{qn,Pn)- 

We then define the kinematical space of states of the quantum theory, Tit, as the space of functions of N real 
variables ip(q) that are square integrable. In this space we define operators Q and P as usual. To construct the 
operators at other time levels (in the Heisenberg Picture) we introduce a linear invertible operator U that we will 
define later and we take 



Q n = U~ 1 Q n -\U = U- n Q U n , P n = U^Pn-iU = U- n P U n . 



(5) 
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When the evolution is determined by a discrete Hamiltonian if, as is the case in the uniform discretizations, the 
evolution operator is given by U = e~ %H l h . Notice that U may also be determined by requiring that the fundamental 
operators satisfy an operatorial version of the evolution equations, 

QnU-UQn+l(Q n ,Pn)=0, A, U ~ UP n+l (Q n , P n ) = 0, (6) 

and this provides a consistency criterion for the construction of U . 

At a classical level if = if and only if the constraints <pi = 0. There exists a natural definition of the physical space 
of the continuum theory that does not require that we refer to the constraint. Since we know that U = exp(—iH/K), 
a necessary condition satisfied by the states of the physical space of the continuum theory, ip £ ifph ys is gi ven by 

Utfj — ip. More precisely the states ip of Hp hys should belong to the dual of a space 4> of functions sufficiently regular 
on Wphys- That is, the states ip £ ^phys satisfy 

ip*U^ipdq = [ ip*ipdq, (7) 



where tp £ $. This condition characterizes the quantum physical space of a constrained continuum theory without 
needing to implement the constraints as quantum operators by using the discretization technique. 

The unitary operators of the discrete theory allow to construct the "projectors" onto the physical space of the 
continuum theory, which is one of the main goals of any quantization procedure based on Dirac's ideas. It should be 
noted that these are really generalized projectors in the sense that they project to a set of functions that belong in 
the dual of a subspace of sufficiently well behaved functions of Hk- All of this is achieved without having to define 
the quantum constraint. To construct the "projectors" one can compute, 

P = lim C M U M . (8) 

M—KX 

If such a limit exists for some Cm such that limAf^ 00 (C , M+i/C'A/) = 1 then UP = P, and we have that UPip = Pip, 

W- £ n k . 

The limit exists in several examples in which H has a continuum spectrum, as we shall see. If the spectrum is discrete 
with eigenvalues et and it contains a vanishing eigenvalue ei then a projector is trivially defined as |e< >< 6j |. A 
constructive procedure leading to a general definition of the projector in terms of the discrete evolution operator U 
valid for any spectrum, continuum or discrete, is given by: 

Int(rAf) - 

P= lim y CnU . (9) 

m->oo Int(rM) -M y ' 

where r is a real number grater than one and Int(rAf) is the integer part of rM . If U has a continuum spectrum 
this definition is a trivial consequence of the previous one. In the case of a discrete spectrum one can check that the 
definition works recalling the definition of the Kronecker delta in terms of a Fourier series. Notice that the definition 
of physical space that Thiemann introduces in the "phoenix project" 0, is equivalent to the choice we make if one 
is considering the Hamiltonian that is quadratic in the constraints. Furthermore, given two states of 7f p hy S , V'ph, 
p h, where ipph — Pip, and 4> p h = P<fi, the physical inner product is defined by < ^ph|0ph >= / dq4>(q)* Pip{q) and a 
physical inner product is determined by the projector constructed from the discrete theory. 

We now illustrate the technique with a rather general example. We consider a generic mechanical system with 
a finite dimensional phase space with one constraint cf> = 0. We will show that the projector constructed with our 
technique reproduces the one constructed with group averaging techniques 0- That is, 

P = lim v&e^ = r — e 1 ^- (10) 

To make contact with the group averaging case we need to assume that is a self-adjoint operator in the kinematical 
phase space with an eigenbasis given by <p\a >= <f)(a)\a > and 1 = J \a >< a\da. The proof of the equivalence is, 



P = P J \a><a\da = J Yon VI^iMe^ 14 ' 2 \a >< a\ (11) 



and noting that limM->oo v 7 AiirMe lMx — 5{x) — -^e l ^ x the proof is complete. For the proof we assumed a 
quadratic form of the Hamiltonian, but can actually be extended to Hamiltonians of the general form we discussed 
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above, computing the integral by steepest descents. The proof can also be extended to systems with N Abelian 
constraints by noting that \im.M^oo('liTTM) N / 2 e~ lMx - x — 5{x) — d/j, N e l ^ x /(2tt n ). This includes important cases 
like gravity in 2 + 1 dimensions, where one can immediately reproduce the results obtained by Perez and Noui [(| via 
group averaging. 

Summarizing, the method of uniform discretizations allows to tackle satisfactorily the three central problems of 
the dynamics of quantum general relativity and provides new avenues for studying numerically classical relativity 
as well. It is based on a set of discretizations generated by Hamiltonians that contain as a particular case the 
quadratic Hamiltonian of Thiemann's "master constraint programme" . The phase space of the continuum theory 
and the physical inner product can be constructed in a straightforward way from the discrete theory and therefore 
provides a generalization of the group averaging extension of the Dirac procedure to systems with structure functions 
in their constraint algebra, like is the case in general relativity. The use of non-quadratic Hamiltonians is possible 
and adds flexibility to the method. The flexibility is crucial, for instance in tackling in an extremely compact and 
straightforward way gravity in 2 + 1 dimensions (we will discuss this example and other issues in a forthcoming paper). 
In other approaches to the dynamics of quantum gravity a major obstacle is the need to define the constraints as 
quantum operators in an unambiguous way. This may be due to the fact that the constraints are only well defined 
on a diffeomorphism invariant space of sets where the constraint algebra is trivial, or, as in the case of the "master 
constraint" since one has only one constraint. This requires an a-posteriori study of each quantization proposal to 
determine if they can reproduce general relativity in some suitable semi-classical limit. This is complex and difficult 
to carry out. One is therefore left with proposals that one is not even sure if they have any connection with the 
theory one desires to quantize. In the discrete approach U = e~ lH / h must implement the discrete classical evolution 
associated to the canonical transformations and therefore one has a correspondence constructive principle as a guide 
that requires that the quantum evolution equations reproduce the Heisenberg equations associated with the classical 
theory Contrary to other methods, one can construct the continuum theory as approximated by a discrete theory in 
the kinematical Hilbert space. This allows the use of operators that can be genuinely used as quantum mechanical 
clocks and is therefore possible to characterize the evolution in a relational way in terms of conditional probabilities 
• The method is therefore devoid of the usual conceptual problems of canonical quantum gravity. 

In conclusion, just like lattice methods did for QCD, the uniform discretizations transform the quantization of 
constrained systems into a computational exercise. The only major hurdle that could stand in the way is that the 
continuum limit may not exist for the case of full general relativity. In that case, this will be a strong indication 
that the theory does not exist. The method allows to incorporate all the benefits of the kinematics of loop quantum 
gravity || and provides an unambiguous avenue to characterize the dynamics and complete the quantum theory. 
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